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ABSTRACT 

We present BVIJK photometry of the black-hole candidate GRO J1655-40 in full 
quiescence. We report a refined orbital period of 2.62191 ± 0.00020 days. The light 
curves are dominated by ellipsoidal variations from the secondary star. We model the 
light curves with an upgraded code which includes a more accurate treatment of limb 
darkening. Previous models containing a large cool disk are ruled out, and indeed our 
data can be fit with a pure ellipsoidal light curve without any disk contribution. In 
general agreement with previous results, we derive a confidence region of the correlated 
quantities inclination and mass ratio, centered on an inclination of 70.2 ± 1.9° and mass 
ratio of 2.6 =b 0.3, resulting in a primary mass Mi = 6.3 ± 0.5 Mq (all 95% confidence). 
The complex limits and errors on these values, and on the possible disk contribution to 
the light curve, are discussed. 

Subject headings: black holes. X-ray sources, binary stars. 

1. Introduction 

The Soft X-ray Transients (SXT) are a sub-class of accreting low-mass binary systems con- 
sisting of a neutron star or black hole primary, accretion disk, and Roche lobe filling secondary 
star. Instabilities in mass accretion cause these sources to undergo episodic outbursts. The X-ray 
emission during outburst indicates the existence of an accreting compact object. During outburst, 
the accretion flow contributes nearly all the optical flux. In quiescence, on the other hand, the X- 
ray luminosity is very low, and the optical light is dominated by the contribution of the secondary 
star. At this time the orbital characteristics of the SXT's may be studied in some detail (van 
Paradijs and McClintock 1995). In particular, the mass function can be measured, which provides 
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a lower limit to the mass of the compact object. Further information about the mass ratio (from 
line broadening) and inclination (from ellipsoidal variability) can be used to determine the mass of 
the compact object (Bailyn ct al. 1998, and references therein) which in many cases is above the 
3Mq upper limit of the mass of a neutron star (Chitre and Hartle 1976). 

While this procedure has been carried out successfully for several systems (see references in 
Bailyn et al. 1998), the modeling of ellipsoidal variations can be complicated by various other 
contributions to the observed light curves, in particular contributions from the accreting material. 
The accretion flow typically contributes both a constant flux offset and random flickering. The 
constant flux contribution of the disk decreases the observed amplitude of the ellipsoidal variations, 
and is thus degenerate with the inclination to first order (there arc potentially observable shape 
differences). In systems for which the disk contribution is significant, not only is the amplitude of 
the variations decreased, but the aperiodic flickering may result in changes in the light curve from 
measurement to measurement, complicating the process of accurately modeling the light curve. 
One proposed method to minimize the disk contribution is the use of infrared light curves. For 
example, Shahbaz, Naylor, & Charles (1994) suggest that the disk contribution to the light curves 
of A0620-00 are minimized in the infrared. They have directly measured the disk contribution in 
the i^-band spcctroscopically, and find it to be less than 27%, smaller than in the optical (Shahbaz, 
Bandyopadhyay, &: Charles 1999a). 

Models used to fit the ellipsoidal variations of SXT's in quiescence have traditionally assumed 
accretion flow in the form of a geometrically thin, optically thick disk (e.g. Orosz and Bailyn 1997, 
hereafter OB97; van der Hooft et al. 1998). The standard thin a-disk model corresponds to one 
allowed solution to the hydrodynamic equations which describe differentially rotating gas flows 
(Shakura and Sunyaev 1973; Frank, King & Raine 1992). The a-disk model has been successfully 
used to model the outburst cycle of cataclysmic variables {e.g. Cannizzo 1993). However, recent 
findings have called into question the application of standard thin disk to model SXT systems. In 
particular, Narayan and collaborators have suggested that the accretion flow in quiescence consists 
of an outer disk and an optically thin advection dominated flow closer to the compact object 
(Narayan, McClintock, & Yi 1996). 

The SXT GRO J1655-40 (Nova Sco 1994) was discovered on July 7 1994 with the Burst 
and Transient Source Experiment aboard the Compton Gamma Ray Observatory (Zhang et al. 
1994). Because GRO J1655-40 has a relatively luminous F star secondary the quiescent light 
curves are dominated by star light to an unusually large extent. As a result, an unusually precise 
fit of the orbital parameters is possible. OB97 report values of Porfe = 2.62157 it 0.00015 days, 
/(M) = 3.24 ± 0.09 Mq, i = 69.50 ± 0.08°, and Q = 2.99 ± 0.08 (where Q is taken to be M1/M2) 
resulting in a primary mass of 7.09 ± 0.22 M©. van der Hooft et al. (1998) derive a refined Porfe = 
2.62168 lb 0.00014 days from photometric and spectroscopic quiescent data taken before the April 
1996 outburst. They use f{M) = 3.16 ±0.15 (Bailyn et al. 1995b), and report a range of secondary 
mass 1.60 < M2 < 3.10 Mq and i of 63.7 ^ i ^ 70.7°, corresponding to a primary mass range 
of 6.29 — 7.60 M0. Shahbaz et al. (1999b) obtained improved spectroscopy, and report f{M) = 
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2.73 ± 0.09 M© and a range of q = M2/M1 of 0.337 - 0.436. They derive a primary mass of 
5.5 — 7.9 Mq. Using Keck spectroscopy, Israelian et al. (1999) found a rotational velocity for the 
secondary star of Vj-ot sini = 93 it 3 km s^^, which implies a mass ratio of Q= 2.48 ±0.16 (assuming 
synchronous rotation and using the Shahbaz et al. mass function). No infrared data have been 
available to date. 

We have obtained BVIJK light curves in order to explore the nature of the quiescent accretion 
flow and to further constrain the orbital parameters of the system. In §2 we present the data and 
reduction techniques. In §3 we discuss the values and limits on the orbital parameters and disk 
contribution. Our results are summarized in §4. 

2. Observations and Reductions 

Photometry of the source was obtained from 1999 July through October with ANDICAM ^ 
on the 1.0 meter Yale telescope at CTIO, operated by the YALO Consortium (Bailyn et al. 1999). 
ANDICAM is capable of observing at an optical and an infrared wavelength simultaneously. The 
source was observed in the Johnson BVI filters with a Loral 2048 x 2048 CCD and the JK infrared 
filters with a Rockwell 1024 x 1024 HgCdTe HAWAII Array. The infrared channel of ANDICAM 
allows for a dither pattern of up to seven positions. An observing sequence consisted of a single 720 
second exposure in B simultaneously with seven dithered 100 second exposures in J, followed by 
a single 360 second exposure in V simultaneously with seven dithered 45 second exposures in K, 
followed by a single 180 second exposure in / simultaneously with 4 dithered 45 second exposures 
in K. We typically performed one observing sequence per night between 1999 July 7 and 1999 
October 30 (HJD 2,451,376 2,451,481). 

The optical images were bias and flat-field corrected using standard IRAF routines. Sky 
corrections for the infrared frames were constructed by median-combining all frames from a given 
night of the same filter. This sky field was subtracted from each image in the set. A flat-fleld was 
constructed by taking the difference between a dome flat taken with dome-lights on, and one taken 
with dome-lights off. Finally, the images from each set were co-added. 

The program DAOPHOT II (Stetson 1987) within IRAF was used to remove near neighbors 
of the program object and all comparison stars. Aperture photometry was then performed on 
the cleaned images. The magnitudes of at least 30 nearby stable stars were combined to form a 
comparison sequence for the differential photometry. Time-series were constructed in each filter 
(see Figure 1). In order to estimate the accuracy of our photometry, we extracted ~ 5 stable 
comparison stars which span the magnitude range of the source. We measured the scatter of each 
comparison star about its mean, and use the average rms of all the comparison stars as our estimate 
for the scatter in the program object photometry. We find the average rms of these comparison 
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stars to be 0.017 mag in the B and V filters, 0.014 mag in the I, 0.023 mag in J and 0.032 mag in 
K. 

Our BVIJK hght curves are presented in Figure 1 phased on our current best period of 
2.62191 lb 0.00020 days (1 a error). They are plotted with a stable comparison star of similar mag- 
nitude in order to demonstrate the scatter in the photometry. The data have been phased following 
the convention of OB97 in which phase corresponds to the closest approach of the secondary 
star to the observer. The light curves are dominated by ellipsoidal variations, demonstrating two 
minima per orbit. We observe a deeper minimum at phase 0.5, which is explained by increased 
gravity darkening near the inner Lagrangian point. The optical light curves are quantitatively in- 
distinguishable from those reported by OB97 and van der Hooft et al. (1998). This result confirms 
that the data obtained in March 1996 was truly in full quiescence, and demonstrates that the source 
has returned to the same state after the outburst which began in April 1996 (Orosz et al. 1997). 
The JK light curves are similar in shape to the optical curves. The increased scatter in the K 
band is consistent with the increased scatter in the comparison stars, and does not require intrinsic 
source variability. 

We find a deviant high point corresponding to HJD ~ 2,451,476.5. In all cases but the V 
band, data from this date lie significantly above the rest of the light curve, possibly due to a flare 
or a hot spot. This is unexpected, given the small contribution of the accretion flow. We removed 
this point for purposes of model fitting. 

In order to find an improved value for the orbital period we combined our data with the 
OB97 quiescent data (Figure 2). We were unable to phase the combined data set on the OB97 
spectroscopic period of 2.62157 ± 0.00015 days. However, OB97 used both quiescent and outburst 
radial velocity data to calculate this period. As pointed out by Shahbaz et al. (1999b), X-ray heating 
of the secondary star can cause the 'effective center' of the secondary to shift, making a sinusoidal 
fit inadequate. The OB97 period is therefore suspect. Using our data set in conjunction with only 
the photometric data set of OB97, we found a preferred photometric period of 2.62191 it 0.00020 
days from the Phase Dispersion Minimization (Stellingwerf 1978) package in IRAF. The Lomb- 
Scargle (Press and Rybicki 1989) and Clean (Roberts, Lehar, and Dreher 1987) algorithms within 
the Period package ^ gave similar results. This period is consistent with the published photometric 
period of 2.62168 ± 0.00014 days of van der Hooft et al. (1998). 



3. Discussion 

Our primary purpose in measuring multiwavelength quiescent light curves of SXT's is to con- 
strain the disk contribution in quiescence. This is not only interesting in itself but provides a crucial 
input into the determination of inclination and therefore to the determination of the mass of the 
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primary as well as other binary system parameters. 

Past work on GRO J1655-40 has consistently found a model with an extended cool disk to 
best fit the data. In particular, the depth of the minimum at phase 0.5 (OB97 convention) in the B 
light curve could not be fit without a grazing eclipse of the star by the disk. Since the size of such 
an eclipse is a very strong function of inclination, very tight limits on i were derived (OB97, van 
der Hooft et al. 1998). At the same time, the small relative contribution of light from the disk in 
the optical wavebands required a very cool outer disk. The current work has improved on previous 
studies in two significant ways. We use an improved modeling code which uses model atmosphere 
tables to compute local intensities, providing significantly improved limb darkening calculations. 
In addition we have obtained infrared data, allowing us to better constrain the contribution of the 
outer parts of the disk, which are cooler than the star. 

We used the recently developed ELC code (Orosz & Hauschildt 2000, hereafter OH2000) 
to model the light curves. This code, which is loosely based on the code described in OB97, 
uses model atmosphere specific intensities rather than specific intensities computed using black 
bodies and a one- or two-parameter limb darkening law. For the present problem we used the 
NextGen grid for cool giants (Hauschildt et al. 1999). The NextGen models are the most 
comprehensive and detailed models available for cool stars. The model atmospheres are computed 
using spherical geometry, rather than the usual plane-parallel approximation. OH2000 show that 
for low gravity stars (log 51 < 3.5) the limb darkening behavior of the extended spherical models 
deviates significantly from a simple linear law. Since the secondary star in GRO J 1655-40 has a 
relatively low gravity (log 3 < 3.5 in general, and formally reaching zero at the Li point itself) the 
consequences of the nonlinear limb darkening are likely to be important. In general, the NextGen 
models predict a lower intensity near the limb of the star. The overall result is that in most cases the 
amplitude of the ellipsoidal light curve is larger when the NextGen intensities are used compared 
to the light curve computed using black body intensities (see OH2000 for a detailed discussion of 
this point). This effect is maximized in the bluest band passes, where a small downward shift in 
temperature makes the biggest decrease in flux. Note that it is precisely the B band lightcurve 
near phase 0.5 which OB97 suggest requires the partial eclipse. 

Our primary conclusion based on modeling the new data with the new code is that the large, 
cool disk and partial eclipse advocated in our previous study (OB97) and that of van der Hooft et 
al. (1998) is no longer viable. Indeed, we find no evidence for any contribution from the accretion 

flow to the optical/IR lightcurves. In Figure 1 we have plotted our data with the large disk model 
(computed with the improved model atmosphere code) which provides the best fit to the OB97 
data. While this model fits our optical data as well as it fits the OB97 data, it fails dramatically 
in the IR, where the large contribution from the cool disk dilutes the ellipsoidal light curves too 
much. In fact, a model with an eclipsing cool outer disk would require an outer disk temperature 
on the order of 100 K as found by van der Hooft et al. (1998). We argue below that such a low 
outer disk temperature is physically implausible, due to the radiation field of the F-star companion. 
Alternatively, all bandpasses are well fit by a model with no disk contribution or eclipses (see Fig. 
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2). The two models (with or without disk) are statistically indistinguishable in terms of their fit 

to the B band at phase 0.5. There is still a small residual discrepancy, but it is similar in size for 
both models, in contrast with the results of OB97 who found that with their cruder modeling of 
limb darkening the partial eclipse model fit significantly better. 

A large cool disk with an outer temperature of 100 K also seems improbable given the proximity 
of an F-star companion. While it is beyond the scope of this paper to solve the detailed equations of 
gas dynamics and radiative transfer to determine the outer disk temperature, it seems reasonable to 
derive an approximate lower limit on this temperature due to the radiation field of the companion 
star. The ratio of flux from the companion star at the point of the disk nearest the Li point to that 
at the surface of the star itself is approximately {Rstar/ Ddiskf where Rstar is tbe effective radius of 
the secondary star calculated using the Eggleton approximation (Eggleton 1983) and D^isk is the 
closest distance between the center of the secondary star and the outer edge of the accretion disk. 
Thus the black-body temperature generated by the star's radiation field when it impacts the outer 
edge of the (presumably optically thick) disk would be 



where Tg^^e is the temperature due to the radiation field of the star at that point and Tgtar is 
the mean surface temperature of the star. Since all the gas in the outermost regions of the disk 
originated on the secondary star and continues to be irradiated by the secondary at least some of 

the time as it orbits the primary, Tg^ige may provide a rough lower limit to Touter- We emphasize 
that the actual value of of Touter depends on the temperature of the gas as it leaves the LI point 
(which may be considerably cooler than the average temperature of the star), the gravitational 
potential energy gained by the gas as it falls toward the disk, and any density change or thermal 
energy transfer which may occur upon the impact of the gas with the accretion disk. 

Figure 3 displays limits on the extent and outer temperature of the quiescent disk. The allowed 
region is bounded at low Tputer by Tedge- An upper limit on Touter can be derived empirically by 
computing the limits on the disk contamination of the ellipsoidal lightcurves. We have done this by 
approximating the accretion flow as a black-body disk with a single temperature equal to that at the 
outer edge. This provides an appropriate limiting case, as any temperature rise toward the center 
of the disk will increase the total flux contribution of the disk resulting in more stringent limits. 
We model the disk flux contribution as a constant dilution factor to the ellipsoidal lightcurves, 
taking into account the system inclination and an assumed outer disk radius. The f-statistic is used 
to compare the resulting lightcurves with our data, and determine whether they are statistically 
different at the one and three a level. Table 1 lists the allowed one and three a disk dilution limits 
for each filter computed from the model fits. The curves shown in Figure 3 are derived from the V 
band data, which provided the tightest constraint on the allowed disk region. Note that the limits 
for the combined light curves are much more severe than for any one filter, both because there are 
many more data points involved and because any contribution from the disk would be maximized 
in one particular bandpass, and thus would have to be much lower in most of the others. Finally, 
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Figure 3 also shows the disk size at which partial eclipses of the star by the disk would occur. 
Note that this is a somewhat porous limit, since it varies with i and Q, but it demonstrates the 
qualitative effect of the lack of evidence for such eclipses. 

Based on the arguments presented above, wc adopted a model which consists only of ellipsoidal 
variations from the Roche-lobe filling secondary. Hence the only free parameters are the mass 
ratio Q and the inclination i. We adopted a spectral type of F6III (Shahbaz et al. 1999b), which 
corresponds to a mean secondary star temperature of T^s = 6336 K (Gray 1992). We had three basic 
binary system observables, namely the photometric light curves in the five filters, the radial velocity 
curve (Shahbaz et al. 1999b), and the measurement of the mean projected rotational velocity of 
the secondary star Vrot sinz (Israelian et al. 1999). All of these observables were obtained when the 
source was in complete X-ray quiescence {L^ < 10^'^Lopt)- Thus the light curves and the velocity 
curve did not suffer any biases caused by X-ray heating. It is well known that tidal distortion of 
the secondary star can cause distortions in the line profiles from which the rotational velocity is 
measured, leading to a systematic error in Vrotsinz (e.g. Kopal 1959; Marsh, Robinson, &; Wood 
1994; Shahbaz 1998; OH2000). However, the measurement of the rotational velocity was obtained 
at the photometric phase (i.e. the closest approach of the F-star), which is the phase when 
the potential distortions to the line profiles are minimized. We computed numerical broadening 
kernels for phase and compared them to the analytic kernel commonly used and found almost no 
systematic biases. We conclude that the Israelian et al. (1999) value of of Vrotsinz represents the 
true mean projected rotational velocity of the F-star. 

Our fitting procedure makes use of all of the binary system observables and their uncertainties. 
We will analyze three different data sets together. First are the photometric light curves reported 
in this paper. There are a total of 227 observations in the five filters (48, 48, 47, 43 and 41 in 

B, V, I, J, and K respectively). The second data set is the radial velocity curve of Shahbaz et al. 
(1999), which has 29 points. Finally, the single measurement of fjotsini = 93 ± 3 kms^^ provided 
by Israelian et al. (1999) is our third data set. We seek a binary system model which best fits all 
three basic observables simultaneously, i.e. we minimize the total x^'- Xtot = Xphoto ~^ XrmRV + Xrot- 

There is some ambiguity in assigning weights to the various observations when fitting several 
quantities simultaneously. In this case, we simply assigned equal weights to each of the three data 
sets. For the data sets to be weighted equally requires that the values of associated with the 
best fit to each one be the same. A three-parameter sinusoid fit to the Shahbaz et al. (1999) radial 
velocities yields Xrv/^ ~ 0.98. We therefore take the quoted uncertainties on each radial velocity 
measurement at face value. We also adopted the la uncertainty in vsini quoted by Israelian et al. 
(1999). In the case of our photometry, however, using the observational uncertainties derived from 
comparison stars resulted in x^/z^ > 1 for our best fits. Therefore we have arbitrarily increased the 
errors used for each point by a factor of 1.27, which yields x^/^ = 1 for the best fits. 

We then defined a grid of points in the Q-i plane where 2.0 < Q < 3.65 in steps of 0.05 and 
66° < i < 74.8° in steps of 0.2°. At each point in this grid the orbital separation a was adjusted 
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to minimize the total x^. The orbital separation is needed to convert the model velocity curve and 
rotational velocity into physical units, and is also needed for the computation of the light curve, 
since the model atmosphere intensities are tabulated as a function of the gravity in physical units. 

Using this procedure, we find a minimum overall x^ot at Q = 2.6 and i = 70.2° (see Table 2). 
This point and contours denoting the one and three a confidence limits are shown in Figure 
4. Combining the allowed regions in Q and i shown in Figure 4 with the Shahbaz et al. (1999b) 
mass function f{M), we find the probability distribution for Mi shown in Figure 5. To obtain 
this curve we divided the Q — i plane into 50 confidence regions corresponding to the probability 
associated with Xtot- Then we selected equal numbers of Q,i pairs in each confidence region, and 
for each Q,i pair chose values of f{M) based on the Shahbaz et al. (1999b) value and error. A 
more detailed discussion of this procedure may be found in Orosz &: Wade (1999). The probability 
density shown in Figure 5 is obtained by calculating an Mi for each point in the {Q, i, f{M)) space, 
and summing the resulting values of Mi weighted by an appropriate likelihood. The likelihood 
Lqif of each combination of Q, i, and /(M) is determined by 

LQif = cxp[-XQi - X/] 

where Xqi is the reduced x^ determined for the relevant values of Q and i as described above, 
and Xf = (/ ~ fm)/f^f, where fm is the value of f{M) measured by Shahbaz et al. (1999b) and 
(7/ is their quoted error. The distribution shown in Figure 5 is then the normalized probability 
distribution created by summing all the values of Mi weighted by their likelihoods. We find a 
preferred Mi value of 6.3 it 0.5 M© (95% confidence). The peak of this distribution is lower than 
the result quoted in OB97 due to the lower value of /(M) determined by Shahbaz et al. (1999b), 
and the distribution is wider than the errors quoted by OB97 because the limits on the inclination 
are less stringent without the requirement for a partial eclipse. 

4. Summary 

We present BVIJK quiescent light curves of the SXT GRO J1655-40. The shape of our 
light curves is consistent with that of previously published results, supporting the assumption that 
the source truly went into quiescence in 1996 before its most recent outburst. We report a refined 
photometric period of 2.62191±0.00020 days. We find that our best fit to the data at all wavelengths 
is provided by a model with no contribution from an accretion flow. For this fit we find a range in 
i and Q shown in Figure 4, which is consistent with, but somewhat wider than, previous findings. 
We also find a probability distribution of Mi shown in Figure 5 where Mi = 6.3 ± O.5M0(95% 
confidence). We have placed limits on the temperature and physical size of any accretion disk 
component. 
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Fig. 1. — Our BVIJK light curves folded on P^rh = 2.62191 days. The open circles represent the 

high point which was removed from the data for analysis purposes. Optical magnitudes were found 
to be consistent with those of OB97. Overplotted arc the best fit OH2000 models to the OB97 
data set. This model required a partial eclipse by a disk with Touter ~ lOOOK. Note the poor fit in 
J and K due to the excess disk contribution in the infrared. 

Fig. 2. — Our light curves (filled circles) are shown with the quiescent light curves of OB97 (filled 
triangles). Overplotted (lines) is the best fit model with no disk contribution. Below is shown a 
comparison stars of similar brightness (filled squares) to the program object and the residuals to 
the fit (filled circles). The scatter in the residuals matches that in the comparison stars, except 
near phase 0.5 in the B where the limb darkening effects are at a maximum. 

Fig. 3. — Upper and lower limits to the size and temperature of a thin disk accretion flow. The thin 
lower line is T^dge which we argue represents a lower limit to the temperature of the outer disk due 
to illumination by the secondary star. The dashed vertical line represents the limit imposed by the 
lack of grazing eclipses due to an infinitely thin disk, at i = 70.2 it 1.9° and Q = 2.6 it 0.3. This 
should be seen as an approximate limit due to its sensitivity to the chosen mass ratio and specific 
shape information. The upper curves represent the la (bold line) and 2>a (dashed line) upper limits 
on spatial extent and temperature of a single temperature black-body disk which is assumed to 
contribute a constant flux dilution to the ellipsoidal light curves. The limits are derived from the 
V band data, which provided the most stringent limits on the allowed disk size and temperature. 

Fig. 4. — Allowed region \n Q — i space. Contours shown are the one and three cr limits imposed 
both by fits to the photometric data and by the constraint imposed by the measured Vrot sin i of 
Israclian ct al. (1999), as described in the text. The region above the dotted line would result in 
X-ray eclipses of the central source, which are not observed. 

Fig. 5. — The probability distribution of Mi, found by weighting the values Mi derived over a 
range of Q, i, and /(M) by an appropriate likelihood, as described in the text. 
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Table 1. Allowed Disk Dilution 



Filter 


la 


3(7 


B 


21% 


28.5% 


V 


11.5% 


17.5% 


I 


19.5% 


25% 


J 


15% 


24.5% 


K 


30% 


46% 
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Table 2. Orbital Parameters of GRO J1655-40 



Parameter 


Result 


Reference 


Orbital Period, photometric (days) 


2.62191 ± 0.0002 


This paper 


K2 (km s-i) 


215.5 ±2.4 


Shahbaz et al. (1999b) 


Vrotsin? 


93.0 ± 3.0 


Israelian et al. (1999) 


Mass Function f{M)/MQ 


2.73 ± 0.09 


Shahbaz et al. (1999b) 


Q 


2.6 ± 0.3 (95% confidence) 


This paper - see Fig. 4 


Inclination i (dcg) 


70.2 ± 1.9 (95% confidence) 


This paper - sec Fig. 4 


Primary Mass Mi/Mq 


6.3 ± 0.5 (95% confidence) 


This paper - see Fig. 5 


Secondary Mass M2/M0 


2.4 ± 0.4 (95% confidence) 


This paper 


Secondary Radius R2/RQ 


5.0 ± 0.3 (95% confidence) 


This paper 
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